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Abstract 

The neutron-neutron scattering length a nn provides a sensitive probe of 
charge-symmetry breaking in the strong interaction. Here we summarize 
our recent efforts to use chiral perturbation theory in order to systematically 
relate a nn to the shape of the neutron spectrum in the reaction ir~d — > nwy. 
In particular we show how the chiral symmetry of QCD relates this pro- 
cess to low-energy electroweak reactions such as pp — > de + v e . This allows 
us to reduce the uncertainty in the extracted a nn (mainly due to short- 
distance physics in the two-nucleon system) by a factor of more than three, 
to < 0.05 fm. We also report first results on the impact that two-nucleon 
mechanisms of chiral order P 4 have on the n~d — ► nwy neutron spectrum. 



1 Introduction 

Quantum chromodynamics (QCD) is almost symmetric under the interchange 
of the up and down quarks. This is called charge symmetry, and is due to the 
fact that the mass difference — m u is much smaller than the QCD mass 
scale A ~ 1 GeV. This symmetry, which is a subgroup of isospin symmetry 
SU(2)v, is well-respected in strong interactions at low energies, but is softly 
broken by quark-mass differences, and also by electromagnetic effects. The 
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relevant dimensionless parameters governing charge-symmetry breaking are 
therefore md ~ rriu and ^ JL , both of which are less than 1%. Although it is 
indeed generally of this small magnitude, charge-symmetry breaking (CSB) 
has many experimentally verifiable effects in hadronic and nuclear physics, 
such as the neutron-proton mass difference, rho-omega mixing, the binding- 
energy difference of mirror nuclei (e.g., 3 He and 3 H), the recently measured 
forward-backward symmetry for np — > dn° [1], and the dd — > an Q reaction [2]. 
For comprehensive reviews of charge symmetry and its breaking see Ref. [3]. 

The difference between the strong-interaction parts of the nn and pp 
scattering lengths (a^ r — a^) is particularly sensitive to CSB. The scattering 
length parameterizes the zero-energy AW scattering phase shift, 5(p), via: 

a^-hm^M, (1) 

p-o p 

where p is the AW relative momentum. Hence the (strong) nn and pp scatter- 
ing lengths would be equal in the limit of exact charge symmetry. Their dif- 
ference is an important quantity for two, somewhat related reasons. Firstly, 
the nucleon-nucleon scattering lengths are unnaturally large compared to the 
pion Compton wavelength. This is indicative of fine tuning in the AW po- 
tential and in consequence the CSB piece of the AW potential has an impact 
on the scattering lengths that is greatly enhanced [3]: 

a Vnn 

Measurements of a pp — a nn therefore provide significant constraints on CSB 
terms in modern phenomenological AW potentials, e.g., AV18 [4]. Secondly, 
when potentials fit to the currently accepted values a nn = —18.59 ± 0.4 fm 
and a pp = — 17.3±0.4 fm [5] are used to make predictions for binding energies 
of mirror nuclei, they very accurately reproduce the experimental binding- 
energy difference of, e.g., the aforementioned 3 H and 3 He [6]. 

Both a nn and a pp must have electromagnetic corrections applied to them 
in order to extract the strong-interaction part. This correction is huge for the 
pp case, but is under good theoretical control. In the nn case the electromag- 
netic correction is due to a magnetic-moment interaction and is ~ —0.3 fm. 

But on the nn side there is an experimental difficulty in obtaining dense 
enough free nucleon targets. There have been some attempts at doing direct 
nn measurements, the most recent one being pursued at the pulsed reactor 
YAGUAR [7]. However, the more promising approaches so far have been 
based on indirect measurements, where final-state neutrons are detected in 
regions of phase space where they have low relative energy and hence observ- 
ables are sensitive to the nn scattering length. 
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Unfortunately, the two most recent measurements employing the nd — > 
nnp reaction for this purpose extract very different a nn values. Thus, a Bonn 
group reported a nn = —16.1 ± 0.4 fm [8], while a group based at TUNL 
claimed a nn = — 18.7 ± 0.7 fm [9], a 4a disagreement. 

However, experiments at different facilities based on the alternative pro- 
cess n~d — > nwy [10,11], have yielded consistent values for many years. Thus 
these results dominate the "accepted" value of a nn quoted above. The scat- 
tering length is extracted by fitting the shape of the spectrum of neutrons 
emitted from the decay of the pionic deuterium atom. The theoretical un- 
certainty in the a nn extracted using extant calculations [12,13] is ~ ±0.3 fm, 
dominated by the uncertainties in the nn wave function at short distances. 

In the present work we revisit these calculations for radiative pion cap- 
ture on deuterium and take advantage of the modern development of effec- 
tive field theory (EFT), in particular chiral perturbation theory (%PT). By 
using an EFT we have consistency between the wave functions and produc- 
tion/capture amplitudes, a recipe to estimate the theoretical error, and we 
can make systematic improvements when necessary. Also, in the case of xPT, 
we gain a close connection to the underlying theory QCD through QCD's chi- 
ral symmetry. In the next section we describe the key elements of our 
calculation of n~d — > nwy, and in Section 3 we present the results already 
obtained in recent publications [14-16], and also provide a first report on 
substantial improvements of these calculations. 

2 Anatomy of the Calculation 

EFTs circumvent the problem of the large QCD coupling constant at low 
energies. Instead one expands amplitudes in the ratio P/A, where P ~ m n 
is a small energy/momentum of the problem and A ~ 1 GeV is the scale 
of chiral- symmetry breaking. This power counting provides a hierarchy of 
quantum-mechanical amplitudes which allows for an systematic organization 
of the calculation. (We count the electron charge e as one power of P/A.) 

In the application of chiral perturbation theory to nuclear processes, 
classes of graphs must be resummed in order to generate the nuclear bound 
states observed in nature. The original proposal for such resummation is due 
to Weinberg [17]. Applying it to the case at hand we see that the amplitude 
for n~d — > nny should be calculated as: 

A={p\6\i/j d ) + {p\T nn G d\ip d ), (3) 

where \ipd) is the deuteron wave function (which is dominated by modes with 
momenta where is applicable), |p) is a plane wave with the observed rel- 
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ative momentum of the two final-state neutrons, p, T nn is the nn rescattering 
amplitude, and G the free nn Green function. 

Meanwhile O is the operator (technically the two-particle irreducible ker- 
nel) governing the transition ir~np — > nwy. Weinberg proposed that O has a 
well-behaved chiral expansion and so can be calculated in %PT. (For a sum- 
mary of the successful application of this idea to electromagnetic processes 
see Ref. [18].) O has one- and two-body pieces, with the one-body part in 
this case beginning at O(P) with the Kroll-Ruderman term for ir~p — > wy. 
Two-body pieces enter at 0(P 3 ). In this work we report on calculations 
obtained from a partial 0(P 4 ) (next-to-next-to-next-leading order = N3LO) 
calculation of O. Our calculation includes all mechanisms at 0(P 3 ) (N2LO), 
but only the dominant 0(P i ) two-body pieces of O. 

2.1 Chirally inspired wave functions 

In order to reach the desired accuracy in the calculation of A, the NN 
wave functions have to be calculated to an order that is consistent with that 
to which O is obtained. Here this means that they must be computed up 
to 0(P 3 ) and thus include the leading- and sub-leading two-pion-exchange 
corrections to the chiral NN potential [19]. The necessary deuteron and nn 
scattering wave functions are derived starting from the asymptotic states, 
given by the asymptotic normalization As and D/S ratio for the deuteron 
and the effective-range expansion for nn scattering. These are integrated in 
from r = oo using the Schrodinger equation with the chiral one- and two- 
pion exchange potentials. Eventually, we reach a region, at r = 1-2 fm, 
where the chiral expansion for the NN potential breaks down. We take the 
simple approach of introducing a cutoff R in this range and assume that the 
potential for r < R is given by a square well whose depth we adjust to enforce 
continuity of the wave function at r — R. This parameterizes and regularizes 
our ignorance of the short-distance NN physics. It is then important to 
ensure that the result is independent of the cutoff R to the order we are 
working, i.e., that the renormalization-group criteria are fulfilled. 

2.2 One-body amplitudes to NNLO 

The chiral one-body amplitudes have been calculated by Fearing et al. [20] 
up to 0(P 3 ), fitting the available 7J9 — > 7i + n and 7i~p — > ^n data via the 
adjustment of xPT low-energy constants (LECs). A preliminary estimate 
of the size of the N3LO one-body amplitude indicates that it has negligible 
influence on 7i~d — > nnj and hence we do not discuss it further here [21]. 
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2.3 Two-body amplitudes to N3LO 
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Figure 1: The two-body diagrams relevant for n~d — > wry. Left: C(P 3 ). 
Right: (9(P 4 ). Only one representative vertex ordering is given for each type 
of diagram. The black discs indicates insertions from C^ N+nNN . 



At 0(P 3 ) (N2L0), there are the three diagrams given to the left in Fig. 1. 
The first is believed to be larger than the second since the pion mass dis- 
appears from its single propagator and the pion can go on-shell, while the 
second has, in addition, also one off-shell propagator. The third diagram 
vanishes in Coulomb gauge. At N3L0 [(9(P 4 )], a slew of diagrams appear, 
given to the right in Fig. 1 and discussed in detail in [16]. The overall result 
for the neutron time-of-flight spectrum when all these graphs except for (/) 
are included is depicted in Fig. 3. Diagram (/) is related to the orthonormal- 
ization of the wave functions. Since this is suppressed by 1/M, we expect 
that effects due to C2-C4, which appear in the two-body currents computed in 
Ref. [16] that are included in our calculation, will be substantially larger. A 
complete calculation of the (9(P 4 ) correction, including orthonormalization, 
is under way [21]. 
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2.4 Constraining unknown short-distance physics 

Fig. 3 shows that the neutron spectrum calculated at 0(P A ) with different 
values of the regulator radius R and a fixed value of the short-distance coef- 
ficient in diagram (to) is significantly different in the final-state-interaction 
(FSI) region. Since a nn is extracted by fitting the shape of the spectrum 
in the FSI region [11] this sensitivity to unconstrained physics of the NN 
system seems to limit the accuracy with which a nn can be obtained from 
it~d — > nny. We now show how to remedy this problem. 

The LO contribution to the matrix element in this region is given by 

r°° f kr\ 

M F si = CJ dru nn (r;p)jo\—-\u d (r), (4) 

where k is the momentum of the outgoing photon, C is a constant that is 
fixed by e, g^, fn, etc., jo is the spherical Bessel function of zeroth order, and 
Ud [u nn (r;p)} is the radial S-wave wave function of the deuteron ( 1 S'o) state. 
But the short-distance part of this matrix element is the same as that of the 
pp fusion matrix element 

POD 

M GT = I dru pp (r)u d (r). (5) 
Jo 

This connection is shown empirically in Fig. 2, revealing a linear relationship 
between the Gamow- Teller matrix element AicT, and the FSI peak height 
in ix~d — > nrry, Tpsi ~ |-Mfsi| 2 - 

This can be understood from the structure of the chiral Lagrangian. It 
contains both one-nucleon and two-nucleon terms linear in the axial field u^. 

C = N^(iv ■ D + g A S -u)N 

- 2d 1 N jf S ■ uNN^N + 2d 2 e abc e KXf ,uV K u x ' a N^ S^r b NN^ S u t c N . . . , (6) 

where 

Uu, = -r a d,n a - e 3ba V,n b r a + f„A, + 0(tt 3 ), (7) 

(A^) is an external vector (axial) field, and the di = O(jjp) are (a 
priori unknown) LECs. Since contains the pion pseudovector coupling, 
as well as a pion-photon coupling and the axial field A^, this is the chiral 
explanation behind the well-known Goldberger-Treiman (GT) relation and 
Kroll-Ruderman (KR) terms. In the two-nucleon sector, the same features 
of Up imply a connection between pion p-wave production, pion photopro- 
duction on the iViV system, and axial currents — two-body analogs of the GT 
and KR. The connection between pion production and electroweak processes 
is currently being investigated by Nakamura 
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Figure 2: Gamow- Teller matrix element plotted against ix~d — > nn7 FSI peak 
height for various values of R and at different orders and wave functions. The 
points correspond to the values of R (in fm) indicated. The straight lines 
are linear fits to the points. The vertical lines show the range in Tfsi after 
renormalization, given the Mqt value indicated by the horizontal line. 



For Gamow- Teller ( 1 S , <-> 3 Si) transitions, the LECs only appear in the 
combination 

d=d 1 + 2(fc + | + ^i + i (8) 

3 3 

where g^d-i = Mf%di and q = Mq [23]. This LEC also appears in p-wave 
pion production in NN collisions, tritium (3 decay, pp fusion, vd scattering, 
fi~d — > nnu^, and the hep reaction. In addition, if the emitted pion couples 
to a third nucleon, this same operator and coefficient enters the leading chiral 
three-nucleon force. These and other implications are discussed further in 
Refs. [15, 16]. The key point in the context of this work is that chiral sym- 
metry and gauge invariance together explain the linear correlation between 
pp fusion and the FSI peak height that is evident in Fig. 2. 

Now, the solar pp fusion process has recently been calculated very accu- 
rately by constraining its unknown short-distance physics from precise cal- 
culations of tritium beta decay [23]. If we adjust the LEC that appears in 
diagram (m) to reproduce this rate for pp — > de + u e we obtain a very precise 
prediction for the FSI peak height in ir~d — > nwy, as shown in Fig. 2. 
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3 Results 




Figure 3: The neutron time-of-flight spectrum for ir~d — > nn7 at different 
cutoffs R and orders as indicated. The thin lines are for the N3LO calculation 
with wave functions calculated with the chiral one-pion exchange potential, 
while the thick lines include N3LO two-body currents and the chiral two-pion 
exchange potential as well. The latter coincide at the quasi-free (QF) peak 
and show a much reduced spread in the FSI peak. 



The result of this renormalization can be seen in Fig. 3. Clearly the 
N3LO contribution reduces the cutoff dependence considerably compared to 
N2LO. The theoretical uncertainty due to unknown short-distance physics 
in the NN system is now negligible in the FSI region. A detailed analysis 
of the other theoretical uncertainties (see Ref. [14]) reveals that the total 
theoretical error in the extracted a nn at N3LO is ±0.3 fm when the entire 
spectrum is fitted and ±0.05 fm if only the FSI peak is fitted. 



4 Conclusions 

Chiral perturbation theory relates the unknown short- distance physics of 
various electroweak two-body observables to pion p-wave production and pion 
photoproduction on two nucleons. (We can also constrain a piece of the chiral 
three-nucleon force from electroweak two-body observables.) This connection 
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makes it possible to calculate n~d — > nn7 to high precision, leading to a small 
theoretical error for the extraction of a nn : Aa^ ory = ±0.05 fm. This reduces 
that error by at least a factor of three compared to previous calculations. 

A future publication [21] will contain a full description of the amplitudes 
and wave functions employed in our N3LO calculation. In that work we will 
also investigate the influence of higher-order electromagnetic corrections in 
the pp wave functions used for pp — > de + u e and whether we are justified in 
neglecting the N3LO one-body contribution. We also provide a full account- 
ing of the 1/M corrections to the two-body operators that are mandated by 
the unitary transformations used to obtain a Hermitian AW potential V. 

In addition we are investigating the possibility to constrain d directly 
from a two-body observable by calculating the \i~ d — > nnu^ capture rate 
in the same framework [24]. This reaction is soon to be measured at the 
Paul Scherrer Institute to 1% precision [25]. It would also be interesting to 
revisit the neutrino-deuteron breakup reactions that are important for the 
SNO detector. Another possible direction would be to complete the circle 
by calculating tritium beta decay using chiral three-nucleon wave functions 
with the r-space regularization we have used in the NN sector. 
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